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Harmonic oscillator in a one-dimensional box 2 

Abstract. We study a harmonic molecule confined to a one-dimensional box with 
impenetrable walls. We explicitly consider the symmetry of the problem for the cases 
of different and equal masses. We propose suitable variational functions and compare 
the approximate energies given by the variation method and perturbation theory with 
accurate numerical ones for a wide range of values of the box length. We analyze the 
limits of small and large box size. 



1. Introduction 

During the last decades there has been great interest in the model of a harmonic 
oscillator confined to boxes of different shapes and sizes [1-23]. Such model has been 
suitable for the study of several physical problems ranging from dynamical friction in 
star clusters [4] to magnetic properties of solids [6] and impurities in quantum dots [23]. 

One of the most widely studied model is given by a particle confined to a box with 
impenetrable walls at — L/2 and L/2 bound by a linear force that produces a parabolic 
potential-energy function V(x) = k(x — xo) 2 /2, where |xq| < L/2. When xq = the 
problem is symmetric and the eigenfunctions are either even or odd; such symmetry is 
broken when Xq ^ 0. Although interesting in itself, this model is rather artificial because 
the cause of the force is not specified. It may, for example, arise from an infinitely heavy 
particle clamped at Xq. In such a case we think that it is more interesting to consider 
that the other particle also moves within the box. 

The purpose of this paper is the discussion of the model of two particles confined to 
a one-dimensional box with impenetrable walls. For simplicity we assume that the force 
between them is linear. In Sec. [2] we introduce the model and discuss some of its general 
mathematical properties. In Sec. [3] we discuss the solutions of the Schrodinger equation 
for small box lengths by means of perturbation theory. In Sec. H]we consider the regime 
of large boxes and propose suitable variational functions. In Sec. |5] we compare the 
approximate energies provided by perturbation theory and the variational method with 
accurate numerical ones. Finally, in Sec. [6] we summarize the main results and draw 
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2. The Model 

As mentioned above, we are interested in a system of two particles of masses rrt\ and 
m 2 confined to a one-dimensional box with impenetrable walls located at x — —L/2 
and x = L/2. If we assume a linear force between the particles then the Hamiltonian 
operator reads 

~ h 2 ( 1 d 2 l d 2 \ k. , 2 

H = -Y{^M + ^dxj) + 2 iXl ~ X2) (1) 
and the boundary conditions are if) — when Xi = ±L/2. It is convenient to convert it 

to a dimensionless form by means of the variable transformation qi = Xi/L that leads 

to: 

where (3 = mi/m 2 , A = kmiL 4 /h 2 and the boundary conditions become ip — if 
qi = ±1/2. Without loss of generality we assume that < (3 < 1. 

The free problem (— oo < < oo) is separable in terms of relative and center-of- 
mass variables 

x = X\ — X2 

X = ~M ( miXl + m 2 x ?) > M = m 1 +m 2 (3) 

respectively, that lead to 

h 2 ( 1 d 2 I d 2 \ . , mi m 2 

F = -T(ii?a5P + mft?J +v ^' n = 4r- (4) 

In this case we can factor the eigenfunctions as 

tpKv(xi,x 2 ) = e tKX (j) v (x) 

-oo < K < oo, v — 0, 1, . . . (5) 

where are the well-known eigenfunctions of the harmonic oscillator and the 

eigenvalues read 

h 2 K 2 Ik ( I 



E "' = ^M + lm{ V+ 2 ■ (6) 
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However, because of the boundary conditions, any eigenfunction is of the form 
ip(xi,X2) = (L 2 /A — x1)(L 2 /A — x 2 ), where $(xi,x 2 ) does not vanish at the 

walls. We clearly appreciate that the separation just outlined is not possible in the 
confined model. 

When P < 1 the transformations that leave the Hamiltonian operator (including 
boundary conditions) invariant are: identity E : ((71,(72) ~ * (91,92) an d inversion 
2 : (91, 92) — ► (—91, —92)- Therefore, the eigenfunctions of Hd are basis for the irreducible 
representations A g and A u of the point group S 2 [24] (also called Ci by other authors). 

On the other hand, when (5 — 1 (equal masses) the problem exhibits the highest 
possible symmetry. The transformations that leave the Hamiltonian operator (including 
boundary conditions) invariant are: identity E : (91,92) — > (91,92), rotation by n 
C 2 '■ (91,92) — ► (92,91), inversion i : (91,92) — > (—91,-92), and reflection in a plane 
perpendicular to the rotation axis a h : (q 1 , q 2 ) — > (—92, — 9i)- In this case the states are 
basis functions for the irreducible representations A g , A u , B g , and B u of the point group 
C 2h [24]. 



3. Small box 



When A < 1 we can apply perturbation theory choosing the unperturbed or reference 
Hamiltonian operator to be H® = H d (X = 0). Its eigenfunctions and eigenvalues are 
given by 

2 cos[(2i - 1)tt9i] cos[(2j - 1)7793] A g 
2 sin(2i7T9i) sin(2j7r9 2 ) A g 

, i,j = 1,2,... 

2 cos[(2i — 1)7791] sin(2j7r9 2 ) A u 
2 sin(2i7T9i) cos[(2j — l)7rg 2 ] A u 

(nl + (3nl), m, n 2 = l,2,.... (7) 

There is no degeneracy when f3 < 1, except for the accidental one that takes place for 
particular values of (5 which we will not discuss in this paper. The energies corrected to 



£ } ,n 2 (9l,92) 



,(0) 
-ni,n 2 



71 

T 



Harmonic oscillator in a one-dimensional box 
first order read 



71 

Y 



A 



2 2 2 o / 2 i 2\ 



(8) 



~ni,na 2 V* ' ' " 127r 2 7l 2 7l| 

When (3=1 the zeroth-order states v^r^ anci Vn^ni ( n i 7^ n 2) are degenerate 
but it is not necessary to resort to perturbation theory for degenerate states in order to 
obtain the first-order energies. We simply take into account that the eigenfunctions of 
adapted to the symmetry of the problem are 



\J%- 5ij {cos[(2i - l)7rgi] cos[(2j - l)7rg 2 ] 
+ cos[(2j - l)7rgi] cos[(2? - l)7rg 2 ]} 

y/2 {cos[(2i - l)7rgi] cos[(2j - l)7rg 2 ] 
— cos[(2j — 1)71-91] cos[(2i — l)7rg 2 ]} 



V2 {cos[(2i - 1)tt9i] sin[2j7rg 2 ] + sin[2j7T9i] cos[(2i - l)7rg 2 ]} A u (9) 



\/2 {cos[(2i — 1)71-91] sin[2j7rg 2 ] — sinpjTrgi] cos[(2i — 1)7T9 2 ]} B u 



^2 — 5ij {sin[2i7rgi] sin[2j7rg 2 ] + sin[2j7rgi] sin[227rg 2 ]} 



\[2 {sin[2i7T9i] sin[2j7rg 2 ] — sin[2j7rgi] sin[2i7rg 2 ]} 
where i,j = 1,2,.... They give us the energies corrected to first order as €$(S) = 



(p\f(S) H d (Pij'iS)^ where S denotes the irreducible representation. Since some of 
these analytical expressions are rather cumbersome for arbitrary quantum numbers, we 
simply show the first sixth energy levels for future reference: 

4" (A, 



A(tt 2 -6) 

7T + 



eF I (4.) = r S- + 



12tt 2 

5tt 2 A(108tt 4 - 405tt 2 - 4096) 



1296tt 4 



e?(B u ) 



5tt 2 A(108tt 4 - 405tt 2 + 4096) 



1296tt 4 



62 ™ ~ 4?r + 247T 2 
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2 A(3tt 2 - 10) 

07T + 



6 

(10) 



*g, -i - • 367r2 

The degeneracy of the approximate energies denoted e^(A g ) and e^(B g ) is broken at 
higher perturbation orders as shown by the numerical results in Sec. El 



4. Large box 

When L — > oo the energy eigenvalues tend to those of the free system (jSj). More precisely, 
we expect that the states with finite quantum numbers n\, n% at L = correlate with 
those with K = when L — ► oo: 



1 + v + 



u = 0,l, 



11' 



ij) Kv {x,X) = < 



(12) 



Besides, we should take into account that the symmetry of a given state is conserved as 
L increases from to oo. 

When /3 < 1 we expect that the states approach 

cos(KX)(j)2v(x) A g 
sm(KX)(p2v+i(x) A g 
sin(KX)(f)2 V (x) A u 
cos(KX)<j) 2v+ i(x) A u 

as L — > oo. 

For the more symmetric case (3 = 1 the states should be 

cos(KX)(t)2 V {x) A g 
sm(KX)4> 2v (x) A u 
cos(KX)4> 2v+ i(x) B u 
sin(KX)(f) 2v+1 (x) By 
Obviously, the perturbation expressions flS} or (fTUl) are unsuitable for this analysis and 
we have to resort to other approaches. 

In order to obtain accurate eigenvalues and eigenf unctions for the present model we 
may resort to the Rayleigh-Ritz variational method and the basis set of eigenfunctions of 
H% given in equations <^ and ([9]). Alternatively, we can also make use of the collocation 



i[) Kv (x,X) = < 



(13) 
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method with the so-called little sine functions (LSF) that proved useful for the treatment 
of coupled anharmonic oscillators [25]. In this paper we choose the latter approach. 

Another way of obtaining approximate eigenvalues and eigenf unctions is provided 
by a straightforward variational method proposed some time ago [26]. The trial functions 
suitable for the present model are of the form 

¥>(fc, ft) = Q - ft 2 ) Q - ft 2 ) /(c, ft, q 2 )e-<^ 2 (14) 

where c = {c\,C2, ■ ■ ■ ,cn} are linear variational parameters, which would give rise to 
the well known Rayleigh-Ritz secular equations, and a is a nonlinear one. Even the 
simplest and crudest variational functions provide reasonable results for all values of A 
as shown in Sec. [5j 

The simplest trial function for the ground state of the model with /3 < 1 is 

«0 = Q - ql) Q - ql) e-<^\ (15) 

Notice that this function is basis for the irreducible representation A g . We calculate 
w(a, A) = (tp\ Hd \(p) I (tp \ f) and obtain A(a) from the variational condition dw/da = 
so that [w(a, A(a)), A(a)] is a suitable parametric representation of the approximate 
energy. In this way we avoid the tedious numerical calculation of a for each given value 
of A and obtain an analytical parametric expression for the energy that we do not show 
here because it is rather cumbersome. We just mention that the parametric expression 
is valid for a > ao where ao is the greatest positive root of A(a) = 0. 

When (3 = 1 we choose the following trial functions for the lowest states of each 
symmetry type 
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Fig. CD shows the ground-state energy for (3 = 1/2 calculated by means of perturbation 
theory, the LSF method and the variational function ( fT5|) for small and moderate values 
of A . Fig. [2] shows the results of the latter two approaches for a wider range of values of 
A. We appreciate the accuracy of the energy provided by the simple variational function 
(I15p for all values of A. The reader will find all the necessary details about the LSF 
collocation method elsewhere [25]. Here we just mention that a grid with N = 60 
was sufficient for the calculations carried out in this paper. Fig. [3] shows that e(A)/\/A 
calculated by the same two methods for (3 = 1/2 approaches \Js/8 as suggested by 
equation (TTTT) . Fig. H] shows the first six eigenvalues e(A) for (3 = 1/2 calculated by 
means of the LSF collocation method. The level order to the left of the crossings is 
ei(Ag) < €i(A u ) < e 2 (AJ < e 2(A g ) < e 3 (A g ) < e 3 (A u ). Notice the crossings between 
states of different symmetry and the avoided crossing between the states 2A U and 3A U . 

Fig. [5] shows the first six eigenvalues for (3 = 1 calculated by means of perturbation 
theory, the LSF method and the variational functions f|T6|) for small values of A. The 
energy order is €i(A g ) < e\(A u ) < e\(B u ) < €2(A g ) < e^{A g ) < £\(B g ) and we appreciate 
the splitting of the energy levels e 3 (A g ) and t\{B g ) that does not take place at first 
order of perturbation theory as discussed in Sec. El Finally, Fig. shows e(A)/v^A for 
sufficiently large values of A. We appreciate that the four simple variational functions 
f fT6|) are remarkably accurate and that e(A)/\/A — > l/y/2 for the first two states of 
symmetry A g and A u and e(A) / \f\ — > 3/\/2 for the next two ones of symmetry B u and 
B g . These results are consistent with equation (|T3|) that suggests that the energies of the 
states with symmetry A and B approach y / 2(2f + 1/2) and v / 2(2f + 3/2), respectively. 
In fact, Fig. [6] shows four particular examples with v = 0. 
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The model discussed in this paper is different from those considered before [1-23] because 
in the present case the linear force is due to the interaction between two particles. 
Although the interaction potential depends on the distance between the particles 
the problem is not separable and should be treated as a two-dimensional eigenvalue 
equation. It is almost separable for a sufficiently small box because the interaction 
potential is negligible in such limit and also for a sufficiently large box where the 
boundary conditions have no effect. It is convenient to take into account the symmetry 
of the problem and classify the states in terms of the irreducible representations because 
it facilitates the discussion of the connection between both regimes. 

The model may be suitable to investigate the effect of pressure on the vibrational 
spectrum of a diatomic molecule and in principle one can calculate the spectral lines by 
means of the Rayleigh-Ritz or the LSF collocation method [25]. 

The simple variational functions developed some time ago [26] and adapted to 

present problem in Sec. H] provide remarkably accurate energies for all values of the box 

size and are, for that reason, most useful to show the connection between both regimes 

and to verify the accuracy of more elaborate numerical calculations. 
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Figure 2. Variational (line) and LSF (points) calculation of the ground-state e(A) for 
0=1/2. 
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Figure 4. First six eigenvalues for /3 = 1/2 calculated by means of the LSF method. 
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Figure 5. First six eigenvalues for {3=1. Circles, dotted line and solid line correspond 
to LSF collocation approach, perturbation theory and variation method, respectively. 
The level order is A g < A u < B u < 2A g < 3A g < B g . 




Figure 6. Variational (solid line) and LSF (symbols) calculation of e(X)/y/\ for j3 = 1. 
The horizontal lines mark the limits l/v2 and 3/v2- 



